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ABSTRACT 

The radar cross section (RCS) of rough rectangular 
plates was studied using Monte Carlo simulations with the 
method of moments (MM) and physical optics (PO), and by 
computing the mean power pattern. For both the MM and PO 
based simulations a plate in the xy plane is represented 
by a collection of small triangles. The roughness of the 
surface is obtained by randomly displacing the z 
coordinates of the triangle nodes. The main focus of the 
study was to look at the effects rough surfaces have on 


the off principal plane regions. 
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I. INTRODUCTION 


Since the invention of radar in World War II many 
applications have been found for its use in both the military 
and civilian sectors. Military radar applications include the 
detection and tracking of ships, aircraft and ground vehicles. 
Modern radar systems make it virtually impossible for ships 
and aircraft to approach targets without being detected at 
great distances. For this single reason, the military has put 
great emphasis on the lowering of the radar cross section 
(RCS) of military vehicles. The word stealth has been used by 
the media to describe specialized military vehicles 
specifically designed with a low RCS. In principle, a stealth 
vehicle can approach a radar without being detected, and thus 
has an increased probability of completing its mission 
objectives relative to a conventional vehicle. 

Several factors contribute to the RCS of a vehicle. The 
shaping of its surfaces and the materials used in its 
construction have the most significant effect. The RCS of a 
rough surface can be significantly different from that of a 
smooth surface. A randomly rough surface reduces specular 
lobes and increases diffuse scattering. Diffuse scattering is 
not desirable because it is isotropic and, even though it has 
a very low level, it cannot be controlled by shaping. Surface 
roughness also affects the behavior of surface waves. 

The sources of target surface imperfections are 
unavoidable. They arise from tolerances in the machining and 
assembly of parts, and in the fabrication of materials 
themselves. For example, the wing of an aircraft typically 
consists of a skin material laid over support structures. The 




skin will have some deviation from the theoretical wing shape. 
Furthermore, if the skin is a composite material, its 
dielectric constant and loss tangent will not be precisely the 
same throughout the material. Thus the wing is a rough, 
inhomogeneous, anisotropic scatterer. 

RCS prediction methods can be divided into two 
categories: 1) approximate methods such as physical optics 
(PO) and 2) rigorous methods such as integral equation 
techniques and the method of moments (MM). In this thesis 
both PO and MM are used to calculate the RCS of plates of 
various sizes, both perfectly flat and rough. For both of 
these methods the plate is represented using triangular 
patches. When MM is applied, the patches are the subdomains 
used in the current expansion. For the PO solution, closed- 
form expressions for scattering from a triangular plate are 
employed. The rough surface is obtained by selecting a random 
number at each node of the plate. The random deviations 
generated in this manner are scaled to produce a random 
surface of a specified variance. The usefulness of PO can be 
accessed by a comparison of the average RCS level with that 
obtained using MM. It is the average RCS level (or diffuse 
scattering) that is of interest in specifying tolerances on 
the manufacturing and assembly processes. 

This thesis investigates the characteristics of 
electromagnetic scattering from surfaces with random roughness 
and its impact on RCS. The first objective is the comparison 
of the approximate and rigorous scattering analysis methods 
using a Monte Carlo simulation. The average RCS levels 
obtained from the Monte Carlo simulations are then compared to 
a closed-form expression for the mean power pattern derived by 
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Ruze (Ruze, 1952). The second objective is to investigate the 
effects of surface roughness on the RCS patterns; in 
particular, how it affects the specular, traveling wave, and 
diffuse RCS components. 

Chapter II is a review of the Rayleigh Criterion and how 
rough surfaces affect specular and diffuse scattering. Also 
discussed is a method for the calculation of the mean power 
pattern using a statistical approach. 

Chapter III reviews PO and MM. The analysis methods are 
presented and the required equations derived. Computer coding 
of the equations in MATLAB is also discussed. 

Chapter IV is an evaluation of the data obtained from the 
computer models, primarily highlighting the differences in RCS 
patterns using the two analysis techniques. A summary 
consisting of the mean power patten difference as related to 
the roughness of the surface is presented. Chapter V 
summarizes the results and highlights the advantages and 
disadvantages of each approach as well as suggests possible 
future work. 
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surface with a plane wave incident from angle 0i in the x-z 
plane. The phase difference (Act)) between the two parallel 
rays scattered into the x-z plane from separate points can be 
calculated by 

A (t>=jt[(Zj-r2)(cos0j+cos0j)+(Xj-Xj)(sin0j-sin0j)] (2.1) 

where the heights of the two separate scattering points are z, 
and Zj. In equation 2.1 the variable k represents the modulus 
of the incident and scattered wave vector (^2n/A) . For 
specular scattering 0 j =02 and the phase difference becomes 

A<|>=2iS:(Zj-Zj)cos0i . (2.2) 

The precise way in which the two scattered rays interact 
with each other depends on the value of AcJ) relative to n. If 
AcJ) « n, then the two scattered waves are nearly in phase and 
thus interfere constructively; if AcJ) = n, then the two 
scattered waves interfere destructively. In the latter case 
there will be no contribution of scattered energy in the 
specular direction. [Ref. 1] The Rayleigh criterion is a 
commonly accepted definition of a rough surface: if A(t> < n 
then the surface is considered smooth, otherwise it is rough. 
For a continuous surface with many scattering points the root 
mean square (RMS) deviation 6 is used in place of Zj-z^. Thus the 
Rayleigh criterion can be expressed as 


R<nM 


(2.3) 







A4) = jt(Xj-x,)(sme,-smej) . 


(2.5) 


As shown in Figure 2, for a smooth plate of very large extent, 
the incident wave is almost completely scattered in the 
specular direction. However, if z, is not equal to the 
surface becomes rough and less energy is scattered in the 
specular direction. Equation 2.2 can be used to find the 
phase difference for the given surface and incident angle. 
When /?„ is small compared to n, then the surface can be 
treated as though it were smooth. If the phase difference is 
no longer negligible, a significant amount of the incident 
energy is scattered diffusely. Finally, if the phase 
variation from point-to-point on the surface becomes large 
then diffuse scattering dominates. A perfectly diffuse 
surface is one that scatters uniformly in all directions for 
every incidence angle. 

B. STATISTICAL APPROACH TO SOLID SURFACE SCATTERING 

Ruze investigated the effects of random surface errors on 
the gain of reflector antennas using a statistical approach. 
[Ref. 2] It is possible to derive a closed-form expression 
for the mean power pattern from the radiation integrals for 
electromagnetic scattering. To derive the expression, some 
simplifying assumptions and approximations are made. One is 
that the surface bumps can be represented locally by a 
Gaussian shape as shown in Figure 3. The width of the 
Gaussion shape is characterized by a correlation interval, c. 
It is the average distance at which the deviation becomes 
uncorrelated. An easy way to explain large correlation 
interval is to imagine a slowly varying deviation as an 
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the mean power pattern can be found in Jenn [Ref. 3] and 
therefore only the highlights are presented here. Along with 
the gaussian approximation just described, a second assumption 
in the derivation of the closed form equation is that, as long 
as the roughness on the surface is sufficiently small, the z 
component of the surface current for a plate lying in the x-y 
plane can be neglected. Only the x and y components contribute 
to the scattered field. For an incident wave polarized 
transverse magnetic to the z-axis (TM,) 


< 2 . 7 ) 

where g=-kT', k=k(ltu+S>v+iw), and r'=itx'andSy'. The quantities u, 
V, and w are the direction cosines of the observation point 
and are defined as: u = sinGsincJ), v = sin9sin(t>, w = cos0, and 
(0,<t)) gives the direction to the observation point. Po 
represents the impedance of free space. 

The scattered field for a unit incident field {Eog=l) is 

£ f ' (2.8) 

where /■„ represents the distance from the local surface origin 
to the observation point. The surface devation as a function 
of position is 6(r'). Now using equation 2.8 in 2.6 
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< P > = < £ •£/ > = //^/ >ds'ds. (2.9) 

The error-free power pattern (i.e. pattern of a prefectly flat 
plate) is 

( 2 . 10 , 

Taking the expected value of the exponential in equation 2.9, 
expanding the integrand in a Taylor Series, then breaking out 
the first term gives 

</.> = ! 'e (2.11) 

where t=t-? and 6^is the variance of the deviations. Switching 
to polar coordinates and extending the limits of the 
integration to infinity for reasons stated by Jenn [Ref. 3], 
equation 2.11 can be further reduced to 

< P ^ ° j| (2.12) 

where A is the area of the plate. 

The RCS, a, can now be found in terms of the mean power 
pattern using the fact that the incident plane wave has unit 
amplitude. For 






Equation 2.10 


be normalized to the area squared 


P 

P =— (2.14) 

and therefore equation 2.13 becomes 


^cos^e 


A 



(2.15) 


The second term in equation 2.15 represents the average or 
diffuse RCS level. Because the variance is present in the 
second term, it shows the importance of maintaining tight 
surface tolerances. 

A similar calculation can be performed for a two 
dimensional array of discrete isotropic scatterers that are 
randomly dislocated from the z=0 plane. The mean power 
pattern consists of two terms similar to equation 2.15. The 
first is the error-free pattern, which depends on the array 
factor of the scatterers. The second term is 


4nA^\ 

CO s^6 (41'co s^0 5 1 

i 

MN j 


where MN is the total number of scatterers making up the 
array. Equation 2.16 represents the average RCS level. It is 
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interesting to note that if the discrete array has a 
rectangular grid spaced dx by dy, its area is MNdyufy. Equating 
the second term of equation 2.15 to equation 2.16 gives 


dxdy=Tic^ . 

Thus an equivalent correlation interval for the discrete array 
can be computed. 
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m. METHOD AND FORMULATIONS 


The equations required for the PO and MM Monte Carlo 
simulations are developed in this chapter. Also introduced is 
a brief description of the computer code used in the 
simulation. 

A. PHYSICAL OPTICS (PO) 

Probably the most commonly used prediction method for the 
study of wave scattering from rough surfaces is the physical 
optics approximation. [Ref. 1] Physical optics provides an 
estimate of the surface current on the illuminated part of the 
target. The shadowed area is assumed to have a surface 
current of zero. One limiting factor that should be mentioned 
is that the use of physical optics will not show the effects 
of traveling waves that normally appear along the edges of 
surfaces. Surfaces that are terminated with edges, as in the 
case of a rectangular plate, have traveling waves reflected 
back towards the source. Thus for monostatic RCS 

calculations, PO can give erroneous results at angles near 
surface grazing. 

The derivation of the scattered fields for physical 
optics begins with the geometrical optics approximation of the 
surface current on the illuminated portion of the body. The 
surface current is 

j ]lfl X H ^ , for illuminated part of the body, (3.1) 

' [ 0, for shadowed part of the body. 

In equation 3.1, H- is the incident magnetic field intensity 
at the surface and fi is the surface normal vector. Now 
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provided that the wave is propagating towards the origin. Now 
substituting the magnetic field intensity into equation 3.1 it 
is possible to approximate the surface current of the 
rectangular plate 

y=-2j5 X (£„g<t)-£ e)-t - . (3.4) 

’lo 

For a flat plate in the x y plane, equation 3.4 becomes 

7=2fx(e£ (3.5) 

[ ’lo J 

where g=x'u+y'v+z‘w. 

Finally, the monostatic scattered field can be determined 
from the surface current 

(3.6) 

Combining equations 3.5 and 3.6 gives 

(3.7) 

23irj ^ 

where A is the total area of the plate and g=h for the 
monostatic case. 

A rough plate can be represented by a collection of 
triangular patches. The edges of adjacent triangles touch, but 
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the normals point in different directions as evident in Figure 
5. Thus the total integral in 3.7 is a sum of integrals, each 
over a single triangle 




(3.8) 


where Nj represents the total number of triangles on the plate 
and /, is the /th integration result. The closed-form 
expressions for /, are derived by Moreira and Prata. [Ref. 4] 
In general, 


iD C,+C. a, C +C. C. ^ 

-2— - F, - —-F.\ (3.9) 

D (D -D) ‘ ^ D (D -D) ^ D D 


where is the area of the triangular cell. Figure 6 shows 
the integration cell geometry used in equation 3.9. For 
plane wave incidence 

C^=A,rAij=0. Co=A,,j. F,=0. 

D=k(B„-BJ. D=k(B„-BJ, Do=kB,3. F^^O. 

Ai_„ and represent the amplitude and phase of the integrand 
at node m of triangle i, respectively, as shown in Figure 6. 
Equation 3.9 can be reduced to 



-^1 


DDJ 


Equation 3.10 has some removable singularities that cause 
problems numerically when I Dpi, |D,| and \ D^-Dp\ are small. 
These singularities are handled as special cases and 
expressions are derived by expanding the integrand in a Taylor 



0.5> 



Figures. Triangular patch model of a rough scattering 
surface. 
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The following speci 
\ <L^ and \D^\iL^ (no 


cases result 



Integration cell geometry in local 









2. If \Dp\<Lc and \D^\<L^ (no restrictions on ID^-Dp] 


‘ fi-0 w-o (n+m+2y. 


If \ Dpi >Lc and |D^|<Lt (no restrictions on |D^-Dp|; 


4. If \Dp\iLt and and |D,-Z?p|<£, 


jD ..0 n\ 


-C^G(n, D^)+C ^-— . 


The function G is defined ; 


G(,n,a)=js"e^'“ds 


I can be evaluated recursively using the relationship 


As noted by Moreira, equation 3.16 is numerically efficient 
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and, accurate results are obtained in most cases by setting 
Lt=.05 and using n, m < 2. However for the rough plate it was 
found that Lt=.001 was required to avoid spikes in the RCS 
pattern. Finally, the RCS of the plate can be calculated 
from 



and by using equation 3.7 and 3.8 to determine the scattered 
field. 

B. COMPUTER SIMULATION USING PHYSICAL OPTICS 

The rough surfaces used for the computer calculations in 
this thesis were created using MBUILD. It is written in 
FORTRAN, and can generate basic surfaces and then approximate 
them with triangular facets. The output is a data file 
consisting of the x, y, and z node locations of each facet. 
The z node is then randomly displaced to generate a rough 
surface of specified standard deviation. For the physical 
optics approximation, this data file is then broken up into 
various MATLAB files consisting of individual x-y-z points, 
edges, vertices, and faces. 

The physical optics approximation calculation was done in 
MATLAB. The code is designed to present the user a three 
dimensional plot of the RCS pattern in direction cosine space. 
Equations 3.10 through 3.14 are functions of the local 
coordinate angles and, therefore, to use these equations, a 
transformation must be performed for each triangle as shown in 
the Figure 7. The transformation involves rotating about the 
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to th 


le plane defined by i and ft, and then rotat; 
ihrough the angle between i and fi in that plane. 


/ Are. 

ft 

< 

/ /i> 



p 

I 




Figure 7. Translation of coordinates for a plate triangle. 


METHOD OF MOMENTS (MM) 

MM is a numerical technique that is very powerful becai 
t can be applied to arbitrary bodies. Basically MM reduc 
n integral equation to a matrix problem. The mati 
[imensions are dependent on the size of the scattering body 
e incident wavelength. For this reason 1 


if th( 





to b( 


lived is the E-field 


sizes. The integral equation 
integral equation (EFIE) 

/ 1, F')-.^[V'-7(F')]V'G(r, F')]^ds ' . (3.19) 

In this chapter, S denotes the surface of integration. 

The MM solution begins by representing the surface current 
by a series with unknown expansion coefficients, I„ : 

■ (3-20) 

J„ are the basis function, and both /„ and J„ can be complex. 
An important consideration in the MM solution is the selection 
of the basis functions. One should chooses a basis function 
that has the behavior expected of the current and minimizes 
computational effort. These functions are typically pulses, 
triangles, sinusoids and 6 functions. 

The next step is to solve for the series coefficients. 
Define a set of weighting functions W„, /7;=1 , 2 , . . . N . It is 
convenient to choose the complex conjugate of the basis 
function being used; this choice is referred to as Galerkin's 
method. [Ref. 3] The EFIE is multiplied by each weighting 
function and integrated over the domain resulting in N 
equations of form 

f /,/^’ 

/|7■a>^7„(F')G(^ F')-—V'-7 „(F')v'G(F, F')| j 
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for m=l,2. ..N. An impedance element can be defined between 
testing function /wand the current basis function n 

(3.22) 

(v-r„(^)jG(F, F') . 

The right side of equation 3.21 is a voltage 

• (3.23) 

Using these definitions, equation 3.21 can be written as 

V„ = Y,I„Z„„ /or m = l,2,../r . (3.24) 

Equation 3.24 can now be presented in matrix form as 

V=Z1 (3.25) 

7 is a vector containing the unlcnown expansion coefficients. 
Z is the impedance matrix and V is the excitation vector. By 
inverting the matrix Z, I can be solved for: 

/=Z‘‘F. (3.26) 

Having evaluated the /„, equation 3.20 is now defined and can 
be used in the radiation integral (equation 3.6). 
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A measurement matrix, R can be defined with elements 

where £, is a unit radiated plane wave. [Ref. 3] 
polarized wave 

Combining equations 3.27 and 3.28 

/{“=/[0 (3.29) 

where the superscript 0 is added to denote the polarization of 
E,. Note that u, v, and w are the direction cosines for the 
observation direction. The same can be done for a $ polarized 
plane wave. The electric field components are 

43tr(, „-i 

and for 

E^=Ae->^±RX . ' 3 . 31 ) 

Referring to Figure 5, the rectangular plate is broken up 
into a collection of triangular faces. The surface current is 
represented by equation 3.20 where ? is a position vector to 
a source point. Figure 8 shows the patch geometry for two 


Using a 


(3.28) 
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triangles with one coiranon edge. The number of basis f 
used is dependent on the number of edges shared wil 
triangles. If a patch has two shared edges, the 
spanned by two basis functions; a patch with three 
edges would be spanned by three basis functions. In F 



Figures. Triangular patch geometry. 


edge n is shared by two triangular patches, T„* and 
areas A„* and A„'. The shared edge has a length of 







f ^»p« - 

—> r '■« 

- 1 2^,* 

"'U„p; _ y.- 

' 

(3.32) 

Some characteristics of these functions are: 

edges that represent body edges, there is no 

of the surface current, 2) the current normi 

edge n is constant and continuous across the 

divergence of exists and has the expressi 

1) for triangular 

normal component 

al to an internal 

edge, and 3) the 

on 

L 

—. r in T* 

L 

- r in r; . 

A„' 

(3.33) 

The MM impedence elements can now be wr 

itten as 

p« 1 

,j (3.34) 

where 



(3.35) 

) rr - 

4j:J(D£ ^ * 

(3.36) 
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where ifThe excitation vector elements are 

(3.37) 

where the incident field is evaluated at the centroid 


(3.38) 

Using the expression for the impedance elements in 3.34, 
the integrals must be evaluated numerically. This should be 
accomplished by transforming to a set of localized coordinates 
on each patch and defining a set of normalized coordinates as 
in the case of PO. Singularity correction of the self term 
[m=n) of the impedance matrix is generally required. 

D. COMPUTER SIMULATION USING MM 

The MM code used to simulate the varying rough surface was 
a code written by Sandia National Laboratories called PATCH. 
[Ref. 5] The output from the code provides the user with the 
currents coefficients and radar cross sections. PATCH uses 
the geometry file generated by MBUILD, which was previously 
discussed. Minor modifications were made to PATCH to 
calculate the RCS pattern of the rough surfaces in direction 
cosine space for comparison with the PO results. 
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rv. SUMMARY OF RESULTS 


A. ROUGH SURFACE PLATE MODELS 

The surfaces considered in all the simulations are 3X by 
3X square plates. This was the largest plate possible given 
computer limitations. A large plate is desirable so that the 
principle plane sidelobes fall off rapidly, thus allowing the 
diffuse RCS of the off principal planes to dominate. Another 
limitation was the computer run time. 

The FORTRAN code MBUILD was used to construct a plate 
comprised of 1800 triangles as shown in Figure 9. The plate 
node coordinates were then altered to produce the rough 
surfaces shown in Figure 10 and Figure 11. Figure 10 has a 
surface deviation uniformly distributed from 0 to O.IX with a 
mean of 0.05X and. Figure 11 has deviation varying from 0 to 
0.2X with a mean of O.IX. The reason for the O.lX and 0.2X 
deviations was to create a significant effect on the RCS 
pattern in the off principal plane regions. The manipulation 
of the z coordinates was done using a FORTRAN program ZCHN 
which generates a uniform random number from 0 to 1 and then 
is scaled based on the desired maximum amplitude deviation. 
Thus three plates are used in the simulation models in order 
to study their effects on RCS patterns. 
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Figured. Flat plate 3 . OX by 3. OX. 
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Figure 11. Rough plate with 0.2X maximum surface deviation. 
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B. METHOD OF MOMENTS RESULTS 

Method of moments is a rigorous solution of the scattering 
equations and therefore includes the effects of traveling 
waves. The software used to calculate the RCS patterns was 
PATCH. Figure 12 shows the RCS pattern for a flat 3X by 3X 
plate. The traveling waves cause the high RCS levels present 
around the edges of the pattern. The main lobe and principal 
plane side lobes are quite noticeable, but the off principal 
plane lobes are not as well defined because of the traveling 
wave. When the surface roughness is increased to O.IX, as in 
Figure 13, the principal plane lobes become less 
distinguishable because diffuse scattering raises the average 
level. Increasing the roughness to 0.2X , as shown in Figure 
14, results in a relatively uniform RCS level everywhere 
except the main lobe. Even the traveling wave effects have 
been reduced by the rough surface. In Figure 15, a cut is 
shown for u=-0.27. It illustrates the effect of the surface 
roughing: the principal plane lobes are decreased and the off 
principal plane lobes are increased. 


35 









50. 




Figure 13. RCS of a rough plate computed using method of 
moments (3X by 3X, variance = 0.817e-3, and 6 „ 3 ,= . 1 X). 


37 


















C. PHYSICAL OPTICS RESULTS 

Using the MATLAB code POCOS, the RCS of 3A by 3X plates 
where plotted in direction cosine space for comparison with 
the MM data. Figure 16 shows a RCS pattern for a flat plate 
using the PO approximation. Because PO does not include 
traveling waves, the RCS levels around the edge of the 
pattern are noticeably lower than those of MM. The 
periodicity of the lobes is apparent when the RCS is plotted 
in direction cosine space. Figures 17 and 18 show the RCS 
patterns for the rough plates (5„„ = 0. IX and 0.2X). The 
effects of random errors are more noticeable in this case than 
they are for MM because of the absence of traveling waves. As 
the surface roughness increases, the off principal plane lobes 
increase, and, if the roughness increases without bound, the 
scattering from the plate would become completely diffuse. A 
comparison of the rough surface patterns are shown in Figure 
19. A cut of the RCS pattern has been taken at u=-0.27. 
Figure 19 shows the effect of the surface roughing. The off 
principal plane levels in some areas increase by as much as 
35dB, while principal plane levels decrease only slightly. 
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Figure 16. RCS of a flat plate computed using physical optics 
(3X by 3X ]. 
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Figure 18. RCS of a rough plate computed using physical optics 
(3A by 3A, variance = 3.271e-3, and 6„3,= .2A). 
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D. COMPARISON OF THE METHODS 

The closed-form approximation, equation 2.15, can be used 
to calculate the average RCS level based on the area, 
variance, and correlation interval of the plate surface. The 
second term in equation 2.15 is plotted in Figures 20 and 21. 
The curves for several correlation intervals are shown. The 
MM and PO RCS curves are averaged over the four quadrants. 
Furthermore values of |u|<0.2 and |v|<0.2 where discarded when 
averaging, because the RCS in these regions is dominated by 
the specular RCS (ie., first term in equation 2.15). Surface 
roughness is not as noticeable in the MM plot because the 
surface wave scattering dominates at wide angles, even off the 
principal planes. The PO patterns illustrate that the diffuse 
scattering does dominate over the specular scattering at wide 
angles (6^40°) . 

The plate nodes are located every O.IX in the x and y 
directions, and the roughness is achieved by randomly choosing 
the z coordinate at each node. Therefore, a first order 
estimate of the correlation interval might be c=0.1X. This 
value gives curves too high relative to PO, which is the best 
indicator of the diffuse level (because it has no traveling 
waves). From the curves shown in the figures, it appears that 
the correlation interval is actually less, perhaps 0.07X. 
This is close to the result obtained equating the second terms 
of equations 2.15 and 2.17. From equation 2.17 


or 


2_ (O.OU)^ 


c=0.06A. . 
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Note that the error terms represent the average levels due to 
surface roughness. The RCS has a Rayleigh probability density 
function and it is expected that ninety-eight percent of the 
sidelobes will be less than this average value plus 6dB. 
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Figure 21. Comparison of the three prediction methods (MM, PO, 
and closed-form) for a 3A by 3X plate and 0.2X. 




V. CONCLUSIONS AND RECOMMENDATIONS 


For both MM and PO, as the surface roughness increased, 
so did the RCS levels in the off principal planes. Because 
MM is rigorous and takes into account traveling waves, it 
represents the most accurate modeling method. It yields 
higher RCS levels in the off principal plane regions than PO. 
However MM has several limitations that should be considered 
before using it. Because MM is so computationally intensive, 
very fast computers with large amounts of memory are required, 
thus limiting the size of object that can be modeled. Another 
factor is the amount of time one is willing to wait for 
results. The plates modeled in this thesis took approximately 
eight days of computer run time per plate on a SparclO with 64 
Megabytes of memory. MM does provide the most accurate 
results but, on the other hand, the PO calculations take 
approximately one day. If an approximate closed-form 
expression could be found for the traveling wave that could be 
added to the PO RCS patterns, the PO approximation would be 
the best choice, yielding results similar to those of MM. 

Ruze's closed-form expression can be used to predict the 
RCS in regions where diffuse scattering dominates. For the 3X 
square plate there is really no angular sector where diffuse 
scattering dominates, and therefore the agreement with MM is 
not good. However, PO does not have a traveling wave lobe and 
therefore diffuse scattering does dominate at wide angles off 
the principal planes. In this case the closed-form results 
agree with PO if the correct correlation interval is chosen. 

Future work should focus on finding an accurate closed- 
form expression for the traveling wave component of RCS that 
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be added directly 


the PO RCS. Thus the 


be expressed approximately by 



That is, the total RCS is a simple summation of specular (PO), 
diffuse (DIF), and traveling wave (TW) terms. Equation 5.1 
assumes that the sum of the fields squared is approximately 
the sum of the powers and this is true if only one of these 
terms dominates in any particular direction. 
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APPENDIX. COMPUTER CODE 


This appendix contains the code developed to calculate RCS 
using the physical optics approximation. Code JPOCOS 
calculates the RCS pattern of a surface of random roughness. 
In present form, the code calculates the scattered RCS levels 
as a function of u and v in increments of .022 or 2/90. This 
creates a data string of 8281 points which can be manipulated 
to create a matrix of 91 by 91 data points. This data is then 
plotted versus u and v giving the three dimensional patterns 
presented in Chapter V. 
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% program POCOS.m 

% latest version 7/31/95 »> Uses geometry data from build <« 

% PO scattering from triangles 

% *** has global to local coordinate transformations *** 

% +♦*♦ vertices must be entered in clockwise order 
clear all 

% Illumination flag: =0 external face only 
iflag=l; 

% Geometry flag for special test cases 
Lt=.0001; % Taylor series region 
Nt=2; % Number of terms in Taylor series 
Rs=.0; % Surface resistivity (normalized) 
wave=l; 
bk:=2* pi/wave; 
rad=pi/180; 

% Incident wave polarization 
Et=l+j*0; %TM-z 
Ep=0+j*0; %TE-z 

% Wave amplitude at all vertices 
Co=l.; 

% Load data arrays generated by build 
load xpts.m 
load ypts.m 
load zpts.m 
nverts=length(xpts); 
load endl.m 
load end2.m 
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nedges=length(end 1); 
load nodel.m 
load node2.m 
load nodeS.m 
nfaces=length(node3 ); 

% Load vind array which gives the three nodes for each triangle 
for i=l:;nfaces 

pts=[nodel(i) node2(i) node3(i)]; 
vind(i,:)=pts; 
end 

ntria=nfaces; 

nvert=nverts; 

x=xpts; 

y=ypts; 

z=zpts; 

% Define position vectors to vertices 
for i=l invert 
r(i,:)=[x(i) y(i) z(i)]; 
end 

% Get edge vectors and normals from edge cross products 
for i=l:ntria 

A=r(vind(i,2), :)-r(vind(i, 1 
B=r(vind(i,3),:)-r(vind(i,l),:); 

C=r(vind(i,3), :)-r(vind(i,2),:); 

N(i,:)=cross(A,B); 

% Edge lengths for triangle "i" 
d(i,l)=norm(A); 
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d(i,2)=norm(B); 

d(i,3)=norm(C); 

ss=.5*sum(d(i,:)); 

Area(i)=sqrt(ss*(ss-d(i,l))*(ss-d(i,2))*(ss-d(i,3))); 

Nn=norm(N(i,:)); 

N(i,:)=N(i,:)/Nn; 
t(i.;)=cross([0 0 l],N(i,:)); 
beta(i)=acos(N(i,3)); 
alpha(i)=atan2(t(i,2),t(i,l)+le-6)-pi/2; 
end 

% Pattern loop 
it=0; 

for u=-l:2/90:l 
for v=-l:2/90:l 
it=it+l; 

Sth(it)=-60; Sph(it)=-60; 
if(u''2+v^2)<=l 

phi=atan2(v,u+1 e-8); 

cp=cos(phi);sp=sin(phi); 

theta=asin(sqrt(u''2+v^2)); 

w=cos(theta); 

st=sin(theta);ct=cos(theta); 

D0=[u vw]; 

uu=ct*cp; vv=ct*sp; ww=-st; 

% Spherical coordinate system radial unit vector 
R=[u V w]; 
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% Incident field in global Cartesian coordinates 
eO(l)=uu*Et-sp*Ep; 
eO(2)=w*Et+cp*Ep; 
eO(3)=ww*Et; 

% Begin loop over triangles 
sumt=0; 
sump=0; 
for m=l:ntria 

% Test to see if front face is illuminated 
ndotk=N(m,:)*R'; 

% if ndotk >=0 

% "local" direction cosines 

ca=cos(alpha(m)); sa=sin(alpha(m)); cb=cos(beta(m)); sb=sin(beta(m)) 
Tl=[ca sa 0; -sa ca 0; 0 0 1]; T2=[cb 0 -sb; 0 1 0; sb 0 cb]; 
D1=T1*D0'; 

D2=T2*D1; 

u2=D2(l); v2=D2(2); w2=D2(3); 

% Find spherical angles in local coordinates 

st2=sqrt(u2''2+v2^2)*sign(w2); % this makes a difference 

ct2=sqrt(l-st2''2); 

phi2=atan2(v2,u2+l e-10); 

th2=acos(ct2); 

cp2=cos(phi2); sp2=sin(phi2); 

% Phase at the three vertices of triangle m; monostatic RCS needs "2" 
Dp=2 ♦bk*((x(vind(m, 1 ))-x(vind(m,3)))*u+... 

(y(vind(m, 1 ))-y(vind(m,3)))*v+... 

(z(vind(m, l))-z(vind(m,3)))*w); 


55 




Dq=2 * bk* ((x(vind(m,2))-x(vind(m,3)))*u+... 
(y(vind(m,2))-y(vind(m,3)))*v+... 
(z(vind(m,2))-z(vind(m,3)))*w); 

Do=2*bk*(x(vind(m,3))*u+y(vind(m,3))*v+z(vind(m,3))*w); 
% Incident field in local Cartesian coordinates (stored in e2) 
el=Tl*eO'; 
e2=T2*el; 

% Incident field in local spherical coordinates 
Et2=e2(l)*ct2*cp2+e2(2)*ct2*sp2-e2(3)*st2; 
Ep2=-e2(l)*sp2+e2(2)*cp2; 

% Reflection coefficients (Rs is normalized to etaO) 
perp=-l/(2*Rs*ct2+l); %local TE polarization 
para=0; %local TM polarization 

if(2*Rs+ct2)~=0 para=-ct2/(2*Rs+ct2); end 
% Surface current components in local Cartesian coordinates 
Jx2=ct2*(-Et2*cp2*para+Ep2*sp2*perp); 
Jy2=ct2*(-Et2*sp2*para-Ep2*cp2*perp); 

% Area integral for general case 
DD=Dq-Dp; 
expDo=exp(j*Do); 
expDp=exp(j*Dp); 
expDq=exp(j*Dq); 

%icase(i,m)=-l; 

% Special case 1 

if abs(Dp) < Lt & abs(Dq) >= Lt 
% icase(i,m)=l; 
sic=0.; 
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for n=0:Nt 

sic=sic+G*Dp)'Wfact(n)*(-Co/(n+l)+expDq*(Co*gfl;n,-Dq))); 

end 

Ic=sic*2*Area(m)*expDo/j/Dq; 

% Special case 2 

elseif abs(Dp) < Lt & abs(Dq) < Lt 
% icase(i,m)=2; 
sic=0.; 
for n=0;Nt 
for nn=0:Nt 

sic=sic+0*Dp)'^*(j*Dq)''nn/fact(nn+n+2)*Co; 

end 

end 

Ic=sic*2*Area(m)*expDo; 

% Special case 3 

elseif abs(Dp) >= Lt & abs(Dq) < Lt 
% icase(i,ni)=3; 
sic=0.; 
for n=0:Nt 

sic=sic+0*Dq)''n/fact(n)*Co*gf(n+l,-Dp)/(n+l); 

end 

Ic=sic*2*Area(m)*expDo*expDp; 

% Special case 4 

elseif abs(Dp) >= Lt & abs(Dq) >= Lt & abs(DD) < Lt 
% icase(i,m)=4; 
sic=0.; 
for n=0;Nt 
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sic=sic+(j*DD)*n/fact(n)*(-Co*gf(n,Dq)+expDq*Co/(n+l)); 

end 

Ic=sic*2*Area(m)*expDo/j/Dq; 

else 

% icase(i,m)=0; 

Ic=2*Area(m)*expDo*(expDp*Co/Dp/DD-expDq*Co/Dq/DD-Co/Dp/Dq); 
end % end of special cases test 
% Scattered field components for triangle m in local coordinates 
Es2(l)=Jx2*Ic; Es2(2)=Jy2*Ic; Es2(3)=0; 

% Transform back to global coordinates, then sum field 
Esl=T2'*Es2'; 

EsO=Tl'*Esl; 

Ets=uu *EsO(l )+w*EsO(2)+ww*EsO(3); 

Eps=-sp*EsO( 1 )+cp*EsO(2); 

% Sum over all triangles to get the total field 
sumt=sumt+Ets; 
sump=sump+Eps; 

% end % end of illumination test 
end % end of triangle loop 
Sth(it)=10*logl0(4*pi*abs(sumt)''2+l.e-10); 
Sph(it)=10*logl0(4*pi*abs(sump)''2+l.e-10); 

end % end of if statement 
end % end of pattern loop V 

end % end of pattern loop U 

Sth(:)=max(Sth(:),-60); 
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Sph(:)=max(Sph(;),-60); 


Sth21=Sth'; 

Sph21=Sph'; 

save Sth21.m Sth21 -ascii 
save Sph21.mSph21 -ascii 
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